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The efficiency of laser-cooling relies on the existence of an almost closed optical-transition cycle 
in the energy spectrum of the considered species. In this respect rare-earth elements exhibit many 
transitions which are likely to induce noticeable leaks from the cooling cycle. In this work, to 
determine whether laser-cooling of singly-ionized erbium Er + is feasible, we have performed accurate 
electronic-structure calculations of energies and spontaneous-emission Einstein coefficients of Er + , 
using a combination of ab initio and least-square-fitting techniques. We identify five weak closed 
transitions suitable for laser-cooling, the broadest of which is in the kilohertz range. For the strongest 
transitions, by simulating the cascade dynamics of spontaneous emission, we show that repumping 
is necessary, and we discuss possible repumping schemes. We expect our detailed study on Er + to 
give a good insight into laser-cooling of neighboring ions like Dy + . 


Rare-earth elements are currently widespread in many 
areas of industry, including telecommunications and elec¬ 
tronics. Over the last years they have also entered 
the field of ultracold matter, for which they present 
very suitable properties M- For example the strong 
magnetic moment of open-4/-shell neutral lanthanide 
atoms, up to 10 Bohr magnetons hb for dysprosium, 
induces anisotropic and long-range dipole-dipole inter¬ 
actions, which make those atoms excellent candidates 
for the production of ultracold dipolar gases @01. In 
particular erbium and dysprosium present bosonic and 
fermionic stable isotopes, for which Bose-Einstein con¬ 
densation and Fermi degeneracy were achieved [11 14]. 


Meanwhile, laser-cooling of trapped ions [15| has al¬ 
lowed for reaching an exceptional control on quantum 
systems, down to the single-particle level 3, 17}. This 


led to the realization of high-precision optical clocks El- 
1251 ]. or of logic gates for quantum-information process¬ 
ing 26-0]. Another noticeable trend in cold matter 
is to merge cold-atom and cold-ion traps, in order to 
study elementary chemical reactions, like charge transfer 
or molecular-ion formation [31B36| . Up to now all laser- 


cooled ions have a similar electronic structure: most often 
one, or possibly two valence electron around a closed- 
shell core, like in alkaline-earth, mercury (Hg + ), ytter¬ 
bium (Yb + ) or indium (In + ) [37|. Laser-cooling of solids 
doped with rare-earth ions was also reported in several 
experiments [0 . 

In this Letter we propose a scheme for laser-cooling 
of open-4/-shell rare-earth ions, taking the example of 
singly-ionized erbium Er + , whose rich electronic struc¬ 
ture yields advantageous properties for ultracold-matter 
physics. For example Er + spectrum is characterized by 
a forest of weak radiative transitions from which emerge 
a few strong transitions 0- 411, adapted to laser-cooling 
and trapping. Moreover the first excited level of Er + , de¬ 
noted G 2 , lies only 440 cm -1 above the ground level Gi 
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FIG. 1. (Color online) (a) Experimentally observed [dj ] en¬ 
ergy levels of Er + sorted by parity, as functions of the total 
electronic angular momentum J. In panels (b) and (c), which 
are zooms on two areas of (a), the levels in solid lines are 
relevant for laser-cooling, while those in dashed fines are not. 


(see Fig. [T] and Table [J); and the G 1 -G 2 transition is al¬ 
lowed both in the electric-quadrupole (E2) and magnetic- 
dipole (Ml) approximations, with widths equal to 0.1 
mHz and 0.9 pHz respectively (see below). The electronic 
structure of the levels Gi and G 2 is so similar, that their 
static dipole polarizabilities only differ by at most 0.2 % 
(see below), which makes the G 1 -G 2 transition weakly 
sensitive to differential Stark shifts 0]. Moreover the 
strong magnetic moment of Er + equal to 8 p-B , opens the 
possibility to observe magnetic-dipole interactions 0]. 

Nevertheless the rich structure of Er + makes it a priori 
difficult to find closed cycles of strong transitions. Ex¬ 
cited ions can decay towards many leaking levels, which 
gives rise to a cascade dynamics that potentially recycles 
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TABLE I. Characteristics of the energy levels i relevant for 
laser-cooling (see Fig. [T] and Ref. 0): total electronic angu¬ 
lar momentum Ji , leading electronic configuration (where the 
xenon core [Xe] has been omitted), leading multiplet, exper¬ 
imental and theoretical energies E°* p and Ej h . The multi- 
plets of levels Gi, G2, Lj, L1-L5 and E4 come from the NIST 
database [45|], while for the other levels, they are extracted 
from our eigenvector calculations. 


label 

i 

Ji 

leading 

configuration 

leading E™ p E\ h 

multiplet (cnR 1 ) (cm -1 ) 

Gi 

13/2 

4/ 1 "( j H 6 )6s(^i /2 ) 

(6, V 2 ) 

0 

14 

g 2 

11/2 

4/ 12 ( 3 7/ 6 )6s( 2 S 1/2 ) 

(6, V 2 ) 

440 

447 

L'i 

15/2 

4/ 11 6s 2 

4 jo 

6825 

6796 

l 2 

17/2 4/ 11 ( 4 // 5/2 )5d6s( 3 D 1 ) 

( 15 / 2 ,1)° 

13028 

12948 

Li 

13/2 

4/ 12 ( 3 R 6 )5d( 2 D 3/2 ) 

(6, 3 / 2 ) 

20728 

20601 

l 2 

13/2 

4/ 12 ( 3 R 6 )5d( 2 D 5/2 ) 

(6, 5 / 2 ) 

22141 

22154 

l 3 

15/2 

4f 12 ( 3 H e )5d( 2 D 3/2 ) 

(6, 3 / 2 ) 

18617 

18655 

l 4 

15/2 

4f 12 CH 6 )5d( 2 D 5/2 ) 

(6, 5 / 2 ) 

21533 

21601 

l 5 

17/2 

4f 12 ( 3 H e )5d( 2 D 5/2 ) 

(6, 5 / 2 ) 

19303 

19302 

Ei 

15/2 4/ 11 ( 4 / 1 ° 5/2 )5d6s( 1 .D 2 ) 

( 15 / 2 ,2)° 

18463 

18412 

e 2 

15/2 4/ 11 ( 4 // 1/2 )5d6s( 3 T> 3 ) 

( 11 / 2 > 3)° 

25712 

25761 

E3 

15/2 

4/ 11 ( 4 /r 5/2 )5d 2 

O 

29473 

29491 

e 4 

15/2 

4/ 12 ( 3 H 6 )6p( 2 P 3 0 /2 ) 

(6, 3 / 2 )° 

29641 

29655 


the ions to the ground level. In order to determine the 
feasibility of Er + laser-cooling, we have thus made ac¬ 
curate quantum-chemical calculations of Einstein-A co¬ 
efficients, characterizing spontaneous emission from ex¬ 
cited level suitable for laser-cooling towards the ground 
level, but also to all possible leaking levels. By insert¬ 
ing the computed A coefficients in a model giving the 
time-dependent fraction of ions in each electronic state, 
we show that the recycling mechanism found in neutral 
erbium [lj is not efficient enough here, and so we iden¬ 
tify possible repumping schemes based on at least two 
auxiliary transitions. 


Our electronic-structure calculations, which were per¬ 
formed with the Racah-Slater method implemented the 
Cowan codes 0, were described in Refs. 44, 47. Iii}. 
They are composed of three steps, (i) Energies and Ein¬ 
stein coefficients are computed ab initio using a Hartree- 
Fock method including relativistic corrections and com¬ 
bined with configuration interaction (HFR+CI). For each 
parity, the calculated energies depend on quantities like 
direct and exchange Coulombic integrals, which are char¬ 
acteristic of each electronic configuration, (ii) These 
quantities are treated as adjustable parameters in order 
to fit the theoretical energies to the experimental ones by 
the least-square method, (iii) Similarly to energies, our 
theoretical A coefficients depend on a restricted num¬ 
ber of radiative parameters that are also adjusted by fit¬ 
ting our A coefficients with available experimental ones 
50j . In the case of Er + , steps (i) and (ii) resulted in 


the detailed interpretation of the energy spectrum given 
in Ref. [44J. For the evaluation of A coefficients (step 


(iii)), the method of [48}, limited up to now to the tran¬ 
sitions involving only the ground level, is now extended 
to transitions involving an arbitrary number of levels in 
both parities. 

The electronic configurations included in our 
model are: [Xe]4/ 12 6s, [Xe]4/ 12 5d, [Xe]4/ n 6s6p and 
[Xe]4/ n 5d6p for even-parity levels, and [Xe]4/ n 6s 2 , 
[Xe]4/ 12 6p, [Xe]4/ n 5d6s, [Xe]4 / 13 and [Xe]4/ 11 5d 2 for 
odd-parity levels, [Xe] denoting the electronic configu¬ 
ration of xenon, omitted in what follows. For example, 
among the computed even-parity levels, 130 were fitted 
to their known experimental counterparts, using 25 
free energetic parameters, giving a 62-cm -1 standard 
deviation. 

After steps (i) and (ii) the level i is described by a 
Cl wave function |i) = ]C P c ip|p), where |p) formally 
represents an electronic configuration. The theoretical 
Einstein coefficients Afj characterizing the probability of 
spontaneous emission from level j to level i can be ex¬ 
panded 



where the configurations p and q are identical, except for 
one electron that hops from subshells n'i' to n£, with 
£' — i = 0, ±1 for electric-dipole (El) transitions. Unlike 
the coefficients a,ij tPq which are specific to each transition, 
the matrix elements of the monoelectronic r-operator are 
common parameters to all transitions. 

The configurations included in our model give rise to 
ten possible r matrix elements: three for [n£-n'£') = 
( 6 s- 6 p) transitions, namely (p, q) = (4/ 12 6s, 4/ 12 6p), 
(4/ n 5d6s,4/ n 5<i6p) and (4/ n 6s 2 ,4/ n 6s6p), three for 


6p-5d, and four for 5d-4/ transitions. In step 
(iii) we fitted the scaling factors (SFs) f m = 
(n£,p\f\n'£',q)/(n£,p\r\n'£',q} liFR between f matrix 
elements and their computed HFR values, with m = 
( n£pn’£'q ), to have the best agreement between exper¬ 
imental and theoretical A coefficients. Among the ten 
scaling factors, those corresponding to 6p-5d transitions 
on one hand, and to 5d-4/ transitions on the other 
hand, were constrained to be equal; they are called /4 
and /*5 respectively. On the contrary, the three scal¬ 
ing factors corresponding to the strong 6 s- 6 p transitions 
can vary independently; they are called / 1 , f 2 and fs 
for (p, q) = (4/ 12 6s,4/ 12 6p), (4/ n 5d6s,4/ n 5d6p) and 
(4/ n 6s 2 ,4/ n 6s6p) respectively. 

We fitted the SFs f m using the experimental Einstein 
coefficients given by Lawler et al. [41]. Due to strong 
differences between A j!) and A/f p , we excluded 17 of the 
418 transitions of Ref. [41|. We checked that none of 
them were involved in the laser-cooling process discussed 
here. To account for the uncertainty of measurements 
we made 100 fits in which all the A^f p coefficients have 
a random value within their uncertainty range. Averag- 
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TABLE II. Einstein-A coefficients characterizing the sponta¬ 
neous emission from laser-cooling levels E 2 -E 4 to the ground 
level Gi and to the main leaking levels L 1 -L 5 . Our theoreti¬ 
cal values are compared with available experimental ones [4~i| . 
The notation (n) stands for x!0 n . 


j (odd) -)• 

e 2 

E3 

e 4 

rf (s- 1 ) 

1.058(6) 

1.533(7) 

1.546(8) 

i 

4 th 

ij 

A e * p 

Ath 

^4 exp 

Ath 

(even) f. 

(s 1 ) 

(s- 1 ) 

(s' 1 ) 

(s- 1 ) 

(s- 1 ) 

Gi 

1.03(6) 

2.83(7) 

1.48(7) 

1.45(8) 

1.49(8) 

Li 

8.58(2) 

5.4(4) 

4.16(4) 

1.4(5) 

1.58(5) 

l 2 

9.93(1) 

- 

6.01(3) 

6.5(4) 

5.03(4) 

l 3 

1.65(3) 

1.9(5) 

1.53(5) 

1.9(5) 

1.57(5) 

l 4 

1.69(3) 

8.6(4) 

3.31(4) 

8.6(4) 

6.15(5) 

Ls 

2.61(4) 

6.6(5) 

2.46(5) 

5.3(6) 

4.25(6) 


ing the best f m obtained for each of the 100 shots, we 
obtain finally the optimal SFs: /1 = 0.887, /2 = 0.798, 
fs = 0.921, /4 = 0.840 and = 0.817, which give a 
standard deviation (see [3, Eq. (15)) a a = 4.18 x 10 6 
s -1 , with 7V par = 5 free SFs, and N\ ev = 401 transitions. 
With those optimal SFs, 40 % of the A coefficient are 
calculated with a precision better than 12 %. A sat¬ 
isfactory comparison of experimental and theoretical A 
coefficients is presented in Table [TT] for a selection of tran¬ 
sitions relevant for laser-cooling. With this optimal set of 
A coefficients we can calculate the scalar static dipole po¬ 
larizability of any level, using Eq. (4) of Ref. [48|, which 
gives, for Gi and G 2 , the same value of 59.4 ag, ao being 
the Bohr radius. 

In order for instance to characterize the G 1 -G 2 transi¬ 
tion, the Cowan codes also allow for calculating Einstein 
coefficients of E2 and Ml transitions. In principle E2 
Einstein coefficients can be determined using the proce¬ 
dure described above, by changing f into f 2 in Eq. ©, 
and taking configurations p and q in the same parity and 
£ — t' = 0, ±1 or ±2. But due to the absence of experi¬ 
mental data we cannot apply step (iii), and we calculate 
the E2 A-coefficients with HFR f 2 matrix elements. 

In order to avoid leaks towards the levels of the lowest 
configuration, including G 2 , we consider cooling transi¬ 
tions with odd-parity upper levels of total angular mo¬ 
mentum Ji = 15/2. We identify 18 possible candidates 
corresponding to transition energy below 30000 cm -1 . 
Among them five levels give rise to closed transitions with 
the ground level, but with small A coefficients. For in¬ 
stance, the energy of the upper level Ei of the strongest 
closed transition is E e * p = 18463 cm -1 , and the corre¬ 
sponding Einstein coefficient is Aq iEi = 1.68 x 10 4 s _1 . 
By contrast, the three strongest transitions, to levels E 2 , 
E 3 and E 4 , which respectively possess 2, 7 and 75 % of 
4/ 12 6p component, and which are therefore characterized 
by A coefficients larger than 10 6 s _1 , will be considered 
in details in what follows. 

Regarding leaking levels, it is remarkable to note that 


the three levels E 2 , E 3 and E 4 behave in a very simi¬ 
lar way. Calculating the branching ratios (BRs) = 
A*//r!; h , where F 4h is the theoretical inverse lifetime of 

L J ' J J _ 

level j given in Table [III we see that levels j = E 2 , E 3 and 
E 4 mainly decay to the ground level i = G 1 , with branch¬ 
ing ratios respectively equal to 97.1 %, 96.7 % and 96.6 
%. Then the main source of leaks is the level L5, with 
r'ij = 2.5 %, 1.7 % and 2.8 %. Levels L 1 -L 4 give BRs 
between 0.1 and 1 %, while levels not mentioned in Table 
ITDgive BRs below 0.001 %. The common feature of levels 
L 1 -L 5 is their dominant 4/ 12 ( 3 i76)5d( 2 T) :) ) char¬ 

acter, where j = 3/2, 5/2 and Ji = 13/2, 15/2 and 17/2, 
which means that the leaks are essentially due to nps/ 2 - 
(n — l)d 3 / 2 , 5/2 transitions, like in alkaline-earth ions. 

The difference between alkaline-earth ions and Er + is 
that electric-dipole transitions are possible from levels 
L 1 -L 5 to so-called secondary leaking levels, including L[ 
and L' 2l which can themselves decay to lower levels, and 
so forth. To determine the efficiency of recycling to the 
ground level, it is therefore necessary to trace out the 
time-dependent population of each level, when all lasers 
are off [l|. To that end, we consider that, at initial time 
t = 0, all the ions are in the excited level E 2 , E 3 or E 4 ; 
then they cascade to lower levels by spontaneous emis¬ 
sion, until they reach a steady state. The fraction A/ of 
ions in level i obeys the system of differential equations 


dXj 

dt 


n 

-r iXi + y ' AijXj, 

j=*+1 


( 2 ) 


where Fj = )Cj=o A?*- Eq. © the index i formally de¬ 
notes the i-th excited level of the ion. Namely, i = 0 cor¬ 
respond to the ground level Gi, i = 1 to the first excited 
level G 2 , and i = 162 to E 4 . In Eq. ©, n is the initially 
populated level, i.e. Xi(t = 0) = 6i n , and at any time 
the ionic population is conserved i.e. ]G ™ =0 Xi(t) = 1 . 
Equation © can be solved analytically, by starting with 
level n, which decays exponentially, X n = exp(—r n t), 
and which acts as a source for n — 1, n — 2, etc.. Finally 
we obtain the general form 


n 

X i (t) = ^2C ij e- r A (3) 

where the time independent coefficients Cij are given by 
the recursion relation 


Cij = E ( 4 ) 

Li i Jfe= l+ i 

defined for i < k < j, and by Ca = Xi(t = 0 ) = Si n . 

Figure [2] shows the time evolution of the fraction of 
ions in some levels relevant in the cascade dynamics, 
for n = 162 (level E 4 ) and considering El transitions. 
Within a fast time scale T/d r; 10 ns, the ions leave the 
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FIG. 2. (Color online) (a) Time evolution of the fraction of 
ions in some levels relevant for laser-cooling (see Fig. [T] and 
Table UJl. in log-log scale (see Eqs. ([3j and <(4j)). Panel (b) is 
a zoom on fractions close to unity, presented in linear scale. 


TABLE III. Percentage (%) of ions in the three most pop¬ 
ulated levels calculated using Eqs. © and © after t = 1 
second, only with El Einstein coefficients, or with El, E2 
and Ml Einstein coefficients (label “E/M”). 



e 2 

El 

E/M 

e 3 

El 

E/M 

E 4 

El 

E/M 

Gi 

99.72 

99.73 

99.23 

99.25 

99.68 

99.71 

G 2 

0.12 

0.12 

0.62 

0.63 

0.12 

0.13 

u 

0.02 

0.03 

0.02 

0.02 

0.03 

0.03 

other 

0.14 

0.12 

0.13 

0.13 

0.17 

0.17 


level E 4 and populate Gi, Li, L 4 and L 5 . An early quasi¬ 
steady state is reached, with populations corresponding 
to the BRs calculated above. After a fraction a mil¬ 
lisecond, the population of those primary leaking levels 
is transferred to secondary ones, some of which, includ¬ 
ing Lf, and to a large extent G 2 , definitively “trap” the 
ions. On the contrary other levels like Lj can empty 
themselves towards the ground level Gi which reaches 
its steady state of 99.7 % after a fraction of second (see 
Fig. 0b)). 

Table [m] gives the population fractions after 1 second, 
when the steady state is reached (see Fig.©, and for the 
three initial levels E 2 , E 3 and E 4 . In view of the slow 
dynamics observed on Fig. [2] we have also included E2 
and Ml Einstein coefficients in our rate equations (0. 
However their influence is not significant after 1 second, 
because when they could play a significant role, e.g. from 
L 1 -L 5 to Gi in the E2 case, El transitions to secondary 
leaking levels are actually faster. Table HTT1 also indicates 
that the recycling process in not sufficient to ensure an 
efficient cooling, since even after 1 second, a noticeable 
fraction of ions, 0.3 %, is trapped in several metastable 
levels, including \J 2 . The ions in Lf, can decay to L[ 
through an E2 transition, but only after 10 6 s. Therefore 
repumping turns out to be necessary. 


Being the main source of leaks from E2, E3 and E4, 
the level L 5 is a natural candidate for repumping, also 
because spontaneous-emission transitions from L 5 would 
drive the ions to odd-parity high-J levels, and so would 
decrease the probability of recycling to the ground level. 
However L 5 cannot be the only “repumped” level, since 
roughly 2 % of the ions would leave the cooling cycle 
at each transition. One possibility would then be to 
use repumping from the other primary leaking levels Li- 
L4. Repumping the five levels L1-L5 would cause less 
than 10~ 5 of loss per transition; but using five repump¬ 
ing lasers seems experimentally unrealistic. 

Another possibility, in addition to level L 5 , would be 
to inject back to the cooling cycle the ions accumulated 
in level G 2 after a few milliseconds (see Fig. 0. How¬ 
ever direct repumping to levels E 2 -E 4 is not possible, 
at least below the electric-octupole approximation. A 
direct transfer from G 2 to Gi, e.g. through a 7 r pulse, 
is also doable, provided that Gi is empty when the 
pulse is applied. The ions could also be repumped from 
G 2 to an auxiliary odd-parity level which preferentially 
decays to the ground level Gi. The level denoted E’ 
(J E - = 13/2, E™ p = 29628 cm" 1 , = 29636 cm" 1 ) 

seems a good candidate, since the corresponding coeffi¬ 
cients are Aq^, = 1.52(7) s -1 , Aq iE , = 1.02(7) s -1 , 
A e - E , = 3.11(6) s - 1 and Aq 2E , = 1.38(6) s -1 . But again 
a small fraction of ions excited in this level would decay 
to undesired leaking levels. 

In this article we have addressed the feasibility of 
singly-ionized erbium (Er + ) laser-cooling, by modelling 
its energy spectrum and Einstein-A coefficients. The 
most promising way is the closed transition to level Ei, 
which is much weaker than the commonly used transi¬ 
tions, but much st rong er than the E2 transition chosen 
to cool Ca + in Ref. [5l| . Regarding the transitions to lev¬ 
els with 4/ 12 6p character, we observe significant leaks to 
levels belonging to the 4/ 12 5g? configuration. Due to the 
passive role of the 4/ electrons in the leaking process, we 
expect our conclusions to be valid for neighboring rare- 
earth ions like dysprosium, holmium or thulium. 
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